Closed form solutions of linear difference 
equations in terms of symmetric products 

Yongjae Cha* 
Johannes Kepler University 
4040 Linz (Austria) 
ycha@risc. jku.at 

January 22, 2013 

Abstract 

In this paper we show how to find a closed form solution for third order dif- 
ference operators in terms of solutions of second order operators. This work is 
an extension of previous results on finding closed form solutions of recurrence 
equations and a counterpart to existing results on differential equations. As moti- 
vation and application for this work, we discuss the problem of proving positivity 
of sequences given merely in terms of their defining recurrence relation. The main 
advantage of the present approach to earlier methods attacking the same problem 
is that our algorithm provides human-readable and verifiable, i.e., certified proofs. 

1 Introduction 

This paper presents an extension of the algorithm solver |[3j |4] |6l that returns closed 
form solutions for second order linear difference equations to third order linear dif- 
ference equations. The solutions that we are looking for are in terms of (finite) sums 
of squares. This is motivated by applying the algorithm for proving inequalities on 
special functions, i.e., on expressions that may be defined in terms of linear difference 
equations with polynomial coefficients. Conjectures about positivity of special func- 
tions inequalities arise in many applications in mathematics and science. Proving them 
usually requires profound knowledge on relations between these special functions. It 
is well known that there exist many algorithms for proving and finding special function 
identities 1301 [8] |20l [iTl . For automated proving of special functions inequalities only 
few approaches exist. Gerhold and Kauers lfT2l [T5l introduced a method that is based 
on Cylindrical Algebraic Decomposition (CAD). This method has been proven to work 
well on many non-trivial examples [13. ,211 . but even though correctness is easy to be 
seen, termination cannot be guaranteed, hence it is not an algorithm in the strict sense. 
A first attempt to clarify the latter issue has been made in 1 16|. One of the features of 
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proofs of special functions identities is that they usually come with a certificate, i.e., 
some easy to check identity that verifies the proof. The CAD-based approach can not 
hope to have a similar certificate in the near future. The method presented here is a 
first step toward human readable proofs of special functions inequalities, although ad- 
mittedly a representation in terms of sums of squares with positive coefficients is not 
expected to exist for any given input. Besides this application, the results presented are 
of independent interest as they provide difference case counterparts to results obtained 
for the differential case [22 , 28 1. 

First we review the available results in the differential case. Let A; be a differential 
field and Ld E k[d], d = d/dx he a linear homogeneous third order differential opera- 
tor. Singer li22 | characterizes when solutions of Ld can be written in terms of solutions 
of a second order operator in k[d]. Van Hoeij ll28l handles the similar problem when 
the coefficients of the second order operator are restricted to k and shows that it will be 
either of the following cases. 

Case 1 Ld is the symmetric square of a second order operator Kd G k[d] 

Case 2 Ld is gauge equivalent to a symmetric square of a second order operator Kd G 

k[d] 

The definitions of symmetric products and gauge equivalence are recalled in sections |23] 
and |2.4| below. The algorithm given in ||28l returns a second order differential operator, 
Kd e k[d/dx], and a gauge transformation in k[d] that sends solutions of the symmet- 
ric squares of Kd to solutions of Ld for Case 2. 

In the differential case, the symmetric square of Ld has order 5 if and only if we 
are in Case 1. In this case, there is a simple formula that gives Kd- Case 2 is equivalent 
to the symmetric square of Ld having order 6 and a first order right-hand side factor 
in k[d] as well as a certain conic of LdClHl Equation 4.2.1]) having a non-zero solution 
in k. Since for k = C{x) this conic is solvable over C(a;), the last condition becomes 
trivial in this case. The algorithm given in [281 in the first step checks the order of the 
symmetric square of Ld to distinguish the cases. 

The difference case behaves differently; here we denote hy D — C{x)[t] the ring 
of linear difference operators, where r denotes the shift operator Example 12 . 1 5 1 shows 
that the cases can not be distinguished according to the order of the symmetric squares 
when the coefficients are in C(a;). To set up a counterpart theorem for difference equa- 
tions, this example shows that we need one more transformation than that in the differ- 
ential case. Furthermore in Case 1, the algorithm for finding the second order operator 
is more complicated than in the differential case. 

Summarizing, the ideas used in the differential case can not be carried over imme- 
diately to the difference case. Furthermore our aim is to have a closed form solution 
of the given input. Hence, if a factorization is found that is not solvable, this fails to 
satisfy our goal. Thus we build on the ideas of the algorithm solver Ii3]|4l|6l. Here 
we say that a function is in closed form if it is a linear combination of elementary 
functions, special functions or hypergeometric functions over C{x). For instance the 
modified Bessel function of the first kind is a closed form solution of the second order 
operator := zr^ — {2x + 2)t + x + z. 
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The algorithm solver returns closed form solutions for second order linear differ- 
ence operators. The main idea of solver is to map the given operator Li to an oper- 
ator L2 of which a solution is known. This transformation is a bijective map, called 
GT- transformation, that sends solutions of Li to solutions of L2. If a closed form solu- 
tion to one of the operators is known, then by means of this transformation the solution 
of the second operator can be constructed. For this purpose a table with second order 
operators including parameters together with characteristic data (local data) has been 
constructed. This local data can be computed for the given operator, the corresponding 
equivalent operator is found by table look-up. Then by comparing parameters of the lo- 
cal data the GT-transformation can be constructed. The characteristic data is described 
in Section[3] To cover the extension described here the table has been extended so that 
we can give closed form solutions of certain third order linear difference operators. 

2 Preliminary 

In this section we introduce notations used in this paper and recall some known facts ID 
I6l l20ll26l about difference operators. Additionally Cases 1 and 2 above are carried over 
to the difference case for algebraic extensions in Theorem 12 . 1 SI below. 

2.1 Ring of difference oprators 

Let D :— C{x) [t] be the ring of linear difference operators with coefficients in C(x), 
where r is the shift operator acting on a: by T(a;) = x + 1. Then Z) is a noncommutative 
ring where 

T • T^^i = r' for i e N, T • / = T(/)r for / e C{x). 

For L = aij{x)T'^ + • • • + ai{x)T + aQ{x) E D with 7^ 0, we say that L has order d 
and write ord(L) = d. If furthermore ^ then L is said to be a normal operator In 
this paper we will assume all operators to be normal. 

The adjoint operator of L is defined by L* ~ X]f=o o.d-i{x + i)t'' . Suppose L = 
M -Nfov some M, TV e D. Then L* = [M ■ N)* = (r'^i • N* ■ t-"^^) ■ M*, where 
di = ord(M) and thus right-hand side factors of L correspond to left-hand side factors 
of L*. We say that a third order operator L is irreducible in D if both L and L* have 
no first order right-hand side factor in D. 

A second order operator K = h2T^ + hiT + E D is, called a. full operator if 
626160 0- Thus, if X is a normal but not full operator, then 61 = 0. 

2.2 Ring of sequences 

Let := {/ I / : N C}. Then an element v G corresponds to a sequence 
V := (f (1), v{2), 1^(3), . . .). C is embedded in as a subring via constant sequences. 
Suppose v\,V2 G C^, then we define 

V1+V2 ■■= {Vi{l) + V2{l),Vi{2) + V2{2), . . .) 

V1V2 := (t;i(l)t;2(l),i'i(2)t;2(2),...). 
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With the above termwise addition and multiplication, C forms a C-algebra. We define 
the action of T on by r(w) := {v{2),v{'S),v{'i), . . .). 

Let S := C^/^ where si ~ S2 if there exists N e N such that, for all i > N, 
si{i) = S2{i)- Then it is easy to verify that s is a unit in S, i.e. s is invertible in S, if 
and only if s € S has only finitely many zeros. If / G C{x), then the image of / in 
S and the action of r on S are well defined. This way we can embed C{x) to S and 
call s e S rational if there exist g{x) G C{x) and N G N such that g{i) = s{i) for all 
i > N. S[t] forms a ring of difference operators and D is embedded in S[t]. 

We say L{v) = for w e S, L = ad{x)T'^ H h ao(a;) £ S[t] if there is no £ N 

such that 

ad{i)v{i + rf) + ad-i{i)v{i + d — 1) + • • • + ao(i)f (i) = for all i > hq. 

Definition 2.1. h G S is called hypergeometric if r — T{h)/h G S \ {0} is rational 
and r is called the certificate ofh. 

If G S is hypergeometric then (t — r) (h) = where r is the certificate of h. We 
define ¥(1) := {u £ S \ L{u) = 0}. 

Theorem 2.2. [20. Theorem 8.2.1] dimc{V{L)) = OTd{L) for a normal difference 
operator L £ D. 

Thus for a normal operator L E D, V{L) forms a C-vector space with a basis 
{v^eS\l<i< ord(L)}. 

2.3 Term equivalence 

Definition 2.3. The symmetric product, M(§)N, of operators M and N £ D is an 
order-minimal and monic operator such that fiv G V{M(§)N) for all fj, G V{M) and 
V G V{N). 

There is a simple formula if one of the operators has order 1. Let L — ad{x)T'^ + 

••• + ••• + ai{x)T + ao{x) G D and r{x) G C(.t). Then 



Thus, (r — a(x))(S)(T — b{x)) = t — a{x)b{x) for any a{x), b{x) G C(a;). 

Suppose L E D and s G S. Then the above formula gives an operator L = 
L(S)(r - s) G S[t] such that V{L) = {hu \ L{u) = 0} where /i G S is a solution of 
T — s. If L(§){t — s) E D then it is easy to see that s is rational. 

Definition 2.4. Li,L2 E D are said to be term equivalent if there exists T = r—r E D 
such that V{L2) = V{Li(^{t — r)), denoted by Li ^2- Such a T is called the term 
transformation from Li to L2. 
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',^{x) Jl {r{x)) for i = 0, . . . , d - 1. 
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If Li and L2 are term equivalent and r — r is the term transformation then V{L2 ) = 
{hv I h e V{t — r),v £ V{Li)}. Suppose Li and a term transformation T are given, 
then L2 can be obtained by ([1]). 

2.4 Gauge equivalence 

Let Li,L2 e Z) be two given operators, where a closed form solution u of Li is known. 
If furthermore an operator G d D can be determined sucht that G{u) is solution of L2, 
then a closed form solution of L2 can be written as a linear combination of shifts of u 
over C{x). Such a transformation G is called a gauge transformation and Li and L2 
are said to be gauge equivalent if such a transformation exists. 

Definition 2.5. Let Li,L2 G D have the same order. G G D is called a gauge 
transformation/rom Li to L2 iffG : V{Li) — > V{L2) is a bijection. 

Note that G is not required to be a normal operator 

Suppose we are given a gauge transformation G where ord(G) > ord(Li). Then 
there exist Q,G G D with ord(G') < ord(Li) such that G — QLi + G. The remainder 
G is also a gauge transformation that acts in the same way as G on V{Li). Hence, 
w.l.o.g., we may assume that ord(G) < ord(ii). 

Let GCRD(L, M) denote the greatest common right divisor of L, M e D. Since 
G is a bijection, any non zero solution u of Li does not map to zero by G. Thus, Li and 
G have no nontrivial common right hand factor, i.e. GCRD(Li,G) = 1. Using the 
extended Euclidean algorithm G,ii G can be determined such that GG+iiii = 1. 
Then GG is the identity on V{Li) and G is an inverse of G that sends V{L2) —?' V{Li) 
bijectively. 

Definition 2.6. Two operators Li and L2 with the same order are called gauge equiv- 
alent if there exists a gauge transformation G : V{Li) — !■ V{L2) and we use the 
notation Li L2. 

Suppose Li ~g L2 where the gauge transformation from Li to L2 is a single term 
operator, c(a;)r" for n < ord(Li). Then r" • Li • t^" is term equivalent to L2 where 
the term transformation from r" • ii • r^" to L2is t — '^(^+^'> ^ 

2.4.1 How to compute tfie gauge transformation 

Suppose Li and L2 are gauge equivalent and G is a gauge transformation from Li to 
L2- Then there is an operator H G D, ord(iJ) < ord(L2) such that H ■ Li = L2 ■ G. 

The algorithm that was used to find the gauge transformation in ||4ll6l l29l works as 
follows: 

1. For given operators Li and L2, we set up the ansatz G Ci{x)T^, 
where the Ci{x) are undetermined coefficients. 

2. Right divide L2 • G by Li and set the remainder to zero. This way we obtain a 
system A of difference equations for the unknown coefficients Ci{x). 
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3. Compute the rational solutions of the system A to determine the values for 
the Ci{x). 

This algorithm was efficient for second order operators, but for operators of order 
three and higher, computing a solution of the system A, we get at Step|2] is very costly. 
Hence in the current implementation we use the new algorithm HOM to compute the 
gauge transformations that give the set of homomorphisms HomD(V"(ii), V^(L2)) 
in D sending V{Li) to V{L2) for any £1,^2 G D. This means in particular that 
we can drop the condition on the orders, ord(Li) — ord(i2)- 

In short, the algorithm HOM works as follows: For L = Si=o'^«(^)^' ^ 
ad{x) = 1, we define the V-adjoint operator ^^^q a(;_i(x+i— 1)t\ Then there 

is a one to one correspondence between Hom(Li, L2) and rational (invariant under the 
difference Galois group) elements of V{Li) (E) V{L2)- We define a space A^(L^, L2) 
that is isomorphic to V{L^) (8) V{L2)- Then rational elements of A4{LY , L2) corre- 
spond bijectively to elements of Hom(Li , £2)- Thus, we compute rational elements of 
Ai{Li, L2)- This is done by working directly with and L2, and we avoid comput- 
ing large operators such as the symmetric product of Li and L2 (whose solution space 
is a homomorphic image of A4{Li, £2)-) 

Note that if Li and L2 are of the same order, then HOM returns exactly the gauge 
transformations. The algorithm HOM is available at http : / / www . rise . jku . a t /people /ycha/ Horn. txt| 
and more details can be found in [|5|. This is joint work of Yongjae Cha and Mark van 
Hoeij . 

2.5 GT-equivalence 

Definition 2.7. Suppose there is a gauge transformation G and a term transformation 
T = T — r{x) such that the composition of G and T, GoT, maps V{Li) to V{L2), i.e. 
G : V{Li(g){T — r{x)) — V{L2). Then Li and L2 are called GT-equivalent, denoted 
by Li L2, and the composition ofG and T is refered to as the GT-transformation 
from Li to L2. 

Suppose there is a map GT which is a multiple composition of gauge transforma- 
tions and term transformations. Then [19. Theorem 3.3.] shows that we can find a 
gauge transformation G and a term transformation T such that GT{V{Li)) = G o 

nv{L,)). 

2.5.1 How to compute the GT- Transformation 

Definition 2.8. Let C be a subfield ofC and r{x) — cpi{xy^ ■ ■ ■pj{xY^ € G{x),for 
some Ci G Z, monic irreducible in Pi{x) € C[x], and let Si £ G equal the sum of the 
roots ofpi{x). 

r{x) is said to be in shift normal form ;/ — deg{pi{x)) < Re(si) ^ 0, for i = 
1, . . . , J. We denote by SNF(r(x)) the shift normalized form ofr{x), which is obtained 
by replacing each Pi{x) by pi{x + ki) for some ki €z Z such that pi{x + ki) is in shift 
normal form. 
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SNF(r(x)) is unique up to the choice of C. In the algorithm given in Section |5] 
we assume C = Q. For L = ad{x)T'^ + • • • + ao{x) £ D, we denote by det(i) the 
determinant of the companion matrix of L, which is {~l)'^ao{x) /ad{x). 

Theorem 2.9. ^ Theorem 2.3.9] Suppose Li L2forLi,L2 G C(a;)[T] where 
C is a subfield of C Then there exists a gauge transformation G G C(.T)[r] from 
Li® {t — r{x)) to L2for some r{x) G C{x) where 

r{xY = SNF(det(L2)/ det(Li)), ord(Li) = d. 

The original statement of the above theorem uses C = C, but the same proof works 
for any subfield C of C. 

Suppose we know that Li L2 for Li,L2 G Q(x)[t] and we want to find 
the GT-transformation. By the above theorem there exists r{x) G Q_{x) such that 
SNF(det(L2)/ det(ii)) = r{x)'^ where d = ord(Li). When dis even Li(S)(T-r(x)) 
or ii(S)(T + r{x)) can be gauge equivalent to L2. Thus the algorithm Hom will return 
a non-empty set for either of the two. Furthermore, when d is odd, Li(§){t — r{x)) is 
gauge equivalent to L2. 

2.6 Symmetric powers of operators 

Given an operator L E D that annihilates a function u, then in order to obtain an 
operator M £ D that annihilates v? we need the symmetric square of L. In this section 
we state some facts about these operators. 

By L®™ we denote the m* symmetric power of L, i.e., we define L®^ — L and 
L®™ = i(s)i®("~i'. K is called a symmetric square root of i if L = K®"^. 

Suppose Kd is a differential operator of order 2 then it is known that the order of 
L®™ is m + 1 [22. Lemma 3.2, (b)]. However the following lemma shows that this is 
not true for difference operators. 

Lemma 2.10. ^29l Lemma 3] Let K = a2{x)T^ + ai{x)T + ao{x) G D. Then 
L ifai{x) ^ then 

K®^ = bsixy + b2{x)T^ + 6i(x)t + boix), 

where 

b3{x) = ai{x)a2{x + l)^a2(x) 
b2{x) = ai{x + l)a2{x){ao{x + 1)02(2:) 
bi{x) — a^ix + \)ax(x)(a\(x + l)ai(a;) 
b^ix) = —a\{x + l)ao(a; + l)ao(a;)^. 

2. ifax(x) = then K®^ = 02 (a;) - aoix)^. 
The formulas above give order-minimal operators for both cases. 



— ai{x + l)ai{x)) 

- ao(x + l)a2(x)) 
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If a full operator K = a2(a;)T^ + ai{x)T + ao(a;) is a symmetric square root of a 
third order operator L, then also K = a2 (x)t'^ — ai {x)t + uq (x) is a symmetric square 
root of L. If li is a solution of K, then (— is a solution of K. We say K and K 
are conjugates if K K where the term transformation is t + 1. 

Solutions of an equation of type 2 are called Liouvillian solutions ll4l [141 l29l . Sup- 
pose ui is a solution of if = a2(a;)r^ — ao(cc) then {ui,U2}, where it2 = {—lYui, 
forms a basis of V{K) and u\ = u|. Also, it is easy to verify that for arbitrary or- 
ders m it holds that K®"^ = 02 (a;)'" + (-l)'"+iao(a;)™ with a similar proof to the 
one of Lemma l2.10l 

Definition 2.11. A second order operator K is called a unity free operator if the solu- 
tion space of K does not admit a basis {vi, W2} such that — V2 for some n e N. 

Let K — (x^ + x) + (2 x + x"^) r + x^ + 3 X + 2. Then a basis of the solution 
space of K is {xw^f^, xwf } where wi and W2 are solutions of + z + 1 in C. Since 
{xWiY = {XW2Y' = x^ for a: e N, if is not a unity free operator 

Lemma 2.12. If K £ D is an irreducible second order operator then it is a unity free 
operator. 

Proof. We prove this by contraposition. Suppose A' G Z) is not a unity free operator 
Then we may assume K is monic and V{K) admits a basis {wi, W2} such that w" = 
for some n e Z. Let rtg S ^>o be the smallest integer that satisfies w"" = Wj", then 
we may assume vi = {u^qY f = {u^gY f for some / £ S where denotes 
noth root of unity and no,a,b are pairwise relatively prime. Thus K — (r^ — (u"^^ + 
^no)'^ + "no'"no)®(^ ~ whcrc r = T{f)/ f. Sincc ii' is an element in D and by 
equation ([T]i, r is in C{x) and this implies K is reducible in D. □ 

Lemma 2.13. Ifv G S,w 7^ satisfies a full second order operator K = b2{x)T'^ + 
61 (x)r + 5o(x) e I? f/ien v is not a zero divisor in S. 

Proof. We will prove that v has only finitely many zeros. Since K{v) ~ there is 
no G N such that 

b2{x)v{x + 2) + bi{x)v{x + 1) + bfi{x)v{x) = (2) 

and 6i(a;) has no poles or roots for all x > n^^i — 1,2. Suppose = v{ni+l) = 

for some ni > no. Then by (|2|i, =0 for all x > ni and this contradicts that v ^ 0. 
Suppose v{n2) — 0, v{n2 + 1) 7^ for some n2 > no. Then again by (|2]), v{x) ^ 
for all a; > n2 + 1. Thus v{x) 7^ for x large enough and hence u is a unit. □ 

Tlieorem 2.14. If L = K®"^ for some irreducible full second order operator K £ D 

then ord(L) = m + 1 

Proof. Let {vi,V2} be a basis of V{K). We will show that then | i = 

0..m} are linearly independent. Suppose there exist Ci in C, not all zero such that 

CmVT + C„_i<'~\'2 + • • • + CovlJ" = 0. 

By Lemma l2.13l V2 is not a unit and since K is irreducible operator, by Lemma l2.12J 

V1/V2 1 for any n e N. Letz vi/v2 G Sand/(y) := Cmy"^+c,n-iy''^^^^ h 
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Co then /(z) = 0, i.e. f{z{x)) = for all x G N. Thus, z{x) £ {c e C | /(c) = 0} 
for all a; e N and vi = ZV2- Suppose z is not a constant sequence. Since K is an 
irreducible full operator in D, it contradicts that vi is a solution of K. Suppose 2; is a 
constant sequence. Then it contradicts that vi are V2 linearly independent. □ 

In the differential case it holds that if the symmetric square of a third order dif- 
ferential operator Ld € 'C{x)[d] has order 5, then Ld ~ K®"^ for some second order 
operator Kd G C(a;) [d] . However, the following example shows that this does not hold 
in the difference case. 

Example 2.15. Let E := 1)t3 + (-28a;3 -42;" - 36 - 84x- 73a;2)T2 + (-69a;- 
77a;3-18-104x2-4a;^-28a;4)T + a;^ + 5x3 + 8a;2 + 4a; G D. ThenoTd{E®^) = 5. 
A solution of E is xlx{l)'^ where Ix{z) denotes the modified Bessel function of the 
first kind. Then the symmetric square roots of E are Ki = + (2 + 2x)^/ x + It — 
\/ x(x + 1) and K2 = — (2 + 2x)^/x + It — x{x + 1), which are not in D. A 
solution ofKi and K2 are ^/xlx{l) and —^/xlx{l), respectively. 

Let B := zt^ - {2x + 2)t + X + z. Then a solution of B is Ix{z). Then E = 
K®^ ~t B®^, but B and Ki are not gauge equivalent in D, i.e, there is no operator 
in D that sends V{B) to V{Ki). Since ^/x is not a solution of any shift operator in 
D, Hll\ Theorem 5.2]and Lemma A.2], Ki is not a symmetric product of B and a 
difference operator in D. 

In the differential case, suppose Ld K®"^, i.e, the solution of Ld can be ob- 
tained by multiplying a hyperexponential term h to the solutions of Let{wi,U2} 
be a basis of the solution space of Kd, then Ld admits a basis of the solution space 
{gu\^ guiU2,gu\}. However, if g is hyperexponential then ^ is also hyperexponen- 
tial. Thus, Ld = k®^ fork & C(a;)[(5] such that = K(i){d-\^) where 5' = -^h. 

However if is a hypergeometric term, ^/h is not guaranteed to be a solution of an op- 
erator in D. 

Definition 2.16. An irreducible operator L is said to be solvable in terms of second 
order in D if it is GT-equivalent to Ki@K2 ■ ■ ■ ®Kd where the Ki 's are irreducible 
and full second order operators in D. 

We need the following Lemma to prove Theorem l2.18l 

Lemma 2.17. Let Ki,K2 G D be full second order operators. If ord{Ki(^K2) = 3 
then we can choose a basis {vi, V2} ofV{Ki), and a basis {wi, W2} ofV{K2), such 
that V1W2 — V2W1. 

Proof. Let {vi,V2} be a basis of V{Li) and {wi,W2} be a basis of V{L2). Since 
ord(ivri (8)1^2) = 3, the C-vector space generated by {wi^i, ^2^1, has 
dimension 3. Then there exists ai, 02, 03 G C, which are not all zero, such that 

V1W2 = aiviwi + a2V2Wi + azV2'W2. 

Suppose ai — a2 — Q and 03 7^ then it contradicts that vi and V2 are linearly 
independent. Likewise, if 02 = 03 = and oi 7^ then it contradicts that wi and 
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W2 are linearly independent. If oi = 03 = and 02 7^ then we have the desired 
form. So, the remaining cases are either only one of the coefficients ai, 02, 03 is zero, 
or all ai, 02, 03 are non-zero. Here, we will prove the case when 02 is the only zero 
coefficient. Let {wi, ■W2} be another basis of V{L2) such that 

( ^ ( ^ "'A ( 

V*2/ \-ai 1 J \W2) 

Then for {wi, W2} we have V1W2 ~ V2W1 as claimed. □ 

Theorem 2.18. Let L be an operator of order 3, irreducible and solvable in terms of 
second order in D. Then L K®'^ for some irreducible full second order operator 
K Cz D and furthermore 

(a) Lr^t K®^ then oid{L®^) = 5. 

(b) if the gauge transformation of L ^gt K®'^ is not a single term operator then 
ord(L®2) = 6. 

Proof. Let L be a third order, irreducible operator that is solvable in terms of second 
order in D. Then by definition (and the restriction of the order), there exist two ir- 
reducible full second order operators Ki,K2 E D such that L ^gt Ki(§)K2- By 
Lemma 12.171 a suitable basis {vi,V2} for Ki, and a suitable basis {wi, W2} for K2 
can be chosen, such that V1W2 = V2W1. Let h = wi/vi — W2/V2, then /i G S and 
wi — hvi,wAw2 ^ hv2- Since {wi, W2} is a basis for an operator in Z), ft, is hyperge- 
ometric and this implies that Ki K2 with term transformation t — r, where r is the 
certificate of the hypergeometric term h. Summarizing, by Lemma 12.101 L ^gt K®"^ 
for some full operator K E D. 

(a) Let {vi , W2} be a basis of V{K) and r — r be the term transformation from K®"^ 
to L and /i be a solution of t — r. Then {/if ^ , hviV2, /1W2} forms a basis of L and thus 
ord(L®2) = 5 by Lemma|231 

(b) Let {vi,V2} be a basis of V{K). Then {G{hvl), G{hviV2), G{hvl)} forms a 
basis ofV{L), where G — C2{x)t'^ + ci{x)t + co{x) g D is a non single term operator 
and /i is a hypergeometric term. Then G{hvl)G{hv2) ^ G{hviV2)'^ and this impUes 
ord(L®2) = 6. □ 

Suppose Ld is a differential operator of order 3 and Ld — K®"^ for some second 
order differential operator Kd- Then it is well known that there exists a formula to 
construct this Kd, see ll22l Lemma 3.4]. The case where only gauge-equivalence holds, 
i.e., Ld ^g K®'^, is more interesting. In [28] third order operators are treated with a 
focus on determining both Kd and a gauge transformation. 

It is possible to implement a similar algorithm for difference equations which re- 
turns the second order operator I\ to which the given L can be reduced to and a gauge 
transformation. However, in the difference case, in order to give a closed form solution 
of K other algorithms need to be applied or a table look-up. Also, even if we are in 
case (|a)l, finding /v is not as simple as in the differential case, in particular if there 
is a parameter included in the input. Morever to distinguish the cases, the symmetric 
square of a third order operator needs to be computed which can become costly if many 
parameters are involved. 
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3 Local data 



The local data that we are using are the valuation growths at finite singularities in C/Z 
and generalized exponents at the point of infinity. This data is invariant under GT 
transformations. In this section, we review the definition and an invariance property 
(Theorem 13. 61 Theorem 13. 14b of local data from ||3] |4] |6] |27]. We omit proofs in this 
paper. 

3.1 Finite singularities 

Valuation growth was first introduced in [|27| and an algorithm to compute it was given 
in the same paper. Let L = adr'^ + • • • + aor'^ e D. After multiplying L from 
the left by a suitable element of C(x), we may assume that the are in C[x] and 
gcd(ao, . . . , ad) = 1. Then q G C is called a problem point of L if g is a root of the 
polynomial ao{x)ad{x — d) and p G C/Z is called a finite singularity of L if L has 
a problem point in p (i.e. p — q + Z, for some problem point q). Let p E C/Z. For 
a, 6 e p C C we say a > 5 iff a — 6 is a positive integer 

Let e be a new indeterminant, i.e., transcendental over C. We define X^iLo '^i (^+ 

e)T' which is obtained by substituting x with x + e in L. The map L h- > defines an 
embedding (as non-commutative rings) of C(a;)[T] in C{x,e)[T]. Hence, if L = MN, 
then = M^N^. 

Definition 3.1. Let a e C(e) and C'[[e\] be the ring of formal power series over C in e. 
The e-valuation We(a) of a at e = is the element of ZUoo defined as follows: ifaj^O 
then Ve{a) is the largest integer m € Z such that a/e™ G C[[e]], and v^{0) = oo. 

We define an ord(L) dimensional C(£)-vector space 

Vp{L,) :={M:p^C(e) | L,{u)^0}. 

Let qi be the smallest root of ao{x)ad{x — d) in p, so qi is the smallest problem point 
in p. Likewise we define qr to be the largest root of ao{x)ad{x — d) in p. If p is not a 
singularity, that is, if gq and ad have no roots in p, then choose two arbitrary elements 
in p and define qi,qr to be those two elements. 

Definition 3.2. For non-zero u e ^^(Lg) and for a,bCzCifb = a + d— 1, where 
d = ord(Le), we define the box-valuation 

v'^{u) — raiT\{v ^{u{m))\m = a, a + 1, . . . , &}. 
Lemma 3.3. With qi , qr chosen as above, we have 

vlZfiu) = vl\Zi{u) for all q E {qi - 1, - 2, - 3, . . .}, 
^iXlli^) - ^tXliu) for all q e {qr + \.,qr + 2,qr + ?>,...} 

We define Ve^i{u) as i;^'Z^({t) which, by Lemma [331 equals the box valuation of 
any box on the left of qi. Likewise we define v^_r[u) as i'^'^^(u). 
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Definition 3.4. Define the valuation growth of non-zero u S Vp{L^) as 

Define the set of valuation growths ofLatp as 

5p(L) = {gp,,{u) I ii e Vp{L,),u^ 0} C Z. 
If L is a first operator operator then gp{L) has only one element. 

Definition 3.5. Let L be a difference operator and p € C/Z, be a finite singularity of 
L. Ifgp(L) has more than one element then p is called an essential singularity. 

The algorithm given in ||27l determines the set 

\ p is an essential singularity of L}. 

Theorem 3.6. [4^. Theorem 1] ff Li and L2 are gauge equivalent then 

max(g^(Li)) = max(g^(L2)) and min(g^(Li)) = min(g^(L2)) 

for all pe C/Z. 

The following lemma is an immediate consequence of Definition [331 

Lemma 3.7. For each p e C /Z, 

max(g^(L®^) = 2max(g^(L)) and min(g^(i®^) = 2min(g^(L)). 

The above theorem only gives invariance under gauge equivalence. To have invari- 
ance under GT-equivalence, we need to define one more set. Suppose Li L2, then 

Li®{t — r{x)) L2 for some r{x) £ C(x). Then 

max(^(L2)) — max(g^(Li)) + d and min(5^(L2)) = min{g^{Li)) + d 

where {d} — 'g^{T — r{x)), d e Z. So dp {L) = max(g^(L)) — min(g^(L)) is invariant 
under GT-equivalence. Thus, for a difference operator L g, we define a set of ordered 
pairs 

ValG := {{p,dp{L)) € C/Z x Z>o | p is an essential singularity of L}. 

3.2 Singularity at infinity 

LetK := C{{t)),x = 1/i be the field of formal Laurent series and = C{{t^/'')) for 
r S N. We define the valuation for a G K as the smallest power of a whose coefficient 
is non-zero and denote it by v{a). This definition can be extended to D = K[t] = 
K[A], where A :— t — 1 denotes the forward difference, by setting 

v{L) = min{w(aj) + i \ L = ao + ■ ■ ■ + CrfA''} 

for any operator L E D. 
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Lemma 3.8. Let L G K[r]. There exists a polynomial P such that for every n G Zwe 
have 

L{t") = + • • • (3) 

where the dots refer to terms of valuation > n + v{L). 

Definition 3.9. IndL(n), the indicial polynomial of L, is the polynomial P{n) in 
LemmaV 



Lemma 9.2 in Q states that if e Z is a root of IndL(n) then there is u G K 
such that L{u) = and v{u) — N . However, there is no one-to-one correspondence 
between solutions of L in K and integer roots of Ind l {n). For this matter, we introduce 
the ring K[l], where I is a solution of t{1) — l=t, see 1 18 1 for existence of I. We extend 
valuation on K to K[/] by: for a = ait'' + • • • G K[l], G C[/], d G Z, and oi ^ 0, 
we let v{a) = d. With this notion we obtain the following theorem which is equivalent 
to ID Theorem 3.2.10] and lHU Lemma 6.1]. 

Tlieorem 3.10. p £ Z is a solution o/Indi(n) if and only if L has a solution u G 
with v(u) — p. 

An immediate consequence of the above theorem is the following corollary. 

Corollary 3.11. If pi and p2 G Z are the solutions of the indicial equations of Li and 
L2, respectively, then pi + p2 is a solution of the indicial equation of Li(§)L2- 

Define the action of r on Kr as: 



r(t-) = i-(l+<)- 



11.1 1.1 -.,2 



= ^^(l-TT-^ + ^-(- + l)^^ (4) 
-4-(- + i)(- + 2)*' + ---)eK.. 

61 r r r 

Since we have defined the action of r on Kr, we can now apply the formula for the 
term symmetric product in ([T) to K,. [t]. Let and be the following subset and 
subgroup, respectively, of K*: 

Er=iaeK; \ a^ce{^ + ^aif^'^),a^eC,ceC*,v e izj, 



00 

Gr^\a£K:\a = l+ ^ a^f /^ a^ G 

z=r4-l 



Now Ej. is a set of representatives for K* /Gr- The composition of the natural maps 
K* — > K* /Gj. — > Er defines a natural map 

Trunc :K;^ Er. 

Let 

Gr = {aeK;\a = l + —t+ ^ Z'', a, G C, m G Z}. 



i— r+1 
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Definition 3.12. Let r G N then for a,b € E^, we say a is r-equivalent to h, a ~r b, 
when a/b Cz Gr- 

Note that a ~r if and only if = br mod -Z with as in the definition of Er, 
Qi — bi for i < r, and c, v matching as well. 

Definition 3.13. Let g E Er for some r E N. We say that g is a generalized exponent 
of L with multiplicity m if and only if zero is a root of Indian) with multiplicity m 
where L = L(§){t — -). We denote by GenExp(L) the set of generalized exponents 
ofL. 

Suppose L = T — r{x) e D then GcnExp(iy) = {Tnmc(r(t))}. 

Tlieorem 3.14. If two operators Li and L2 are gauge equivalent then for each gi € 
GenExp(Li) there is a g2 E GenExp(i2) such that (72 is equivalent to gi. 

This theorem has been proven first in [6 1. An alternative proof can be found in JSj. 
Tlieorem 3.15. Suppose L,L' E D then 

GenExp(L(S)L') = {Trunc(5g') | g E GenExp(L),g' E GenExp(L')} 

Proof L(s)L'®{t - j^) ^ L(s){t - |)(DL'(D(t - ^7) and since is a solution of 
L(S)(t— i) and L(S)(t— Ois also a solution of the indicial equation of i(S)i'(S)(T— 
T7ini^)byLemma|3il] □ 

Likewise for the valuation growth, we need to define one more set to have invari- 
ance for GT-equivalence. Suppose Li(g)(r — r(x)) L2 for some r{x) E C{x). 
Then 

GenExp(_L2) = {5r3 I 9 G GenExp(Li), {g^} = GenExp(r — r{x))}. 
Thus we define the following set, 

Gquo(L) := {Trnnc{gi/gj) \ g^ ^ gj,gi,gj G GcnExp(i)} 
and then Gquo(ii) = Gquo(L2) if Li ^gt L2. 

4 Table of base equations 

In |l3]|6l, we have formed a table of base equations of order 2, call it TB, as follows; 

• collect equations with known solution from HI |2l . 

• for any closed form expression that shows up often in the literature, generate a 
base equation with existing algorithms lISl flTl . 
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For the algorithm given in Section|5] we have computed symmetric squares of each base 
equation in TB yielding an entry in TB2 of a base equations of third orders. Moreover 
we have generated further base equations as follows: 

Suppose u{x) is a solution of an operator L = J2i=o '^i{x)T^- Then u{x/m) is a 
solution of the operator 

d 
1=0 

As input for our algorithm we accept only operators of order three and the above equa- 
tion may be of higher order One way of obtaining the base equation for u(x/m) in this 
case is using L(„j) when it is a multiple of an operator A/ G D for which M(u(x/m)) = 
0. Since as constructed above is not guaranteed to be the minimial order operator 
we compute Hom(L(^), i(„j)). If the algorithm HOM returns the identity map this 
means that is in fact order- minimal. These cases are neglected and we use L(„j) 
as a base equation only if HOM returns a non-trivial map. 

For instance for the squared hypergeometric function in the table below, 2^1 
an annihilating operator L(2) can be obtained starting from an operator L(i) annihilat- 

using Q. Then the order-minimality of L(2) is checked with 



-x/2+a, x/2+b. 

5 ^ 



2 



2 

—x+a. x+b 



ing 2F1 

the algorithm HOM. In this case HOM returns a non-identity map and hence we save 
L(2) in the table. 

If ct" f e GenExp(L), then zJYrmic{g^f) £ GenExp(L(„j-)), where Zc is a 
root of z"^ — c and e GenExp(T™ — ( — )"). Thus, we can detect whether an 
input operator may have a solution u{x/m) if a base equation for u{x) is in our table. 
However, it is more efficient to compute the base equation for small values of m. 

4.1 Example of base equations 

Here we list a small part of the table which is needed in Section |5] and [T] In the fol- 
lowing table they are listed under (a) a solution (b) the corresponding Gquo, and (c) the 
ValG. The full table can be found at |http : / / www . rise . jku . at /people/ ycha/TB2 . txt| 



1. (a) 2F1 



'x/2+a, x/2+b. 

1 



(b) i^-l,-{2z-l±2^/7^^^f ,±{2z-l±2^/7^^^)^ 

(c) {(-26, 2), (2a, 2), (2a - 2c, 2), (2c - 25, 2)} 
(a) Pxiz)^ (Legendre polynomials squared) 



(b) { 



-l + 2z2±2^/^i^Tl^, (-l + 2z2±2y^i^T^)-\^l±|^!||^^ 



(c) {(0,4)} 

(a) Hx{z)'^ (Hermite polynomials squared) 

(b) {-1 ± V'^2l2T + z^T^, 1 ± 2 V'^2l2T - 4 z^T^} 

(c) {(0,2)} 
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5 Algorithm 

The basic structure of the algorithm is the same that was given in |[3]. Here we use 
an extended table of base equations and a more efficient algorithm for computing the 
gauge transformation, as mentioned in Section l2.4.1l 
Suppose L is the input operator with local data 



for a, &, c e C. By comparing the corresponding data in TB2, we can find that local 
data of L matches with the data of (|2| in Section|4T| Let Ligd be the operator of which 
PxizY is a solution. To compute the parameter z, we compare a with each entry of 
Gquo(Ligd) and compute the set of candidates of possible values for z which is. 



Substituting z by each of the values of the above set, a set of equations cdd2 is ob- 
tained. It remains to cheek for each of the equations in cdd2 whether there is a GT- 
transformation to L and if so then we return the closed form solution by applying the 
GT- transformation io Px{zY. 

Algorithm solver! 

Input: A third order normal operator L/ G Q[2;, r]. 

Output: Either at least one closed form solution of L in the form of cq{x)u{x)'^ + 
ci{x)u{x + 1)^ + C2{x)u{x + 2)^ where Ci{x) are hypergeometric terms and is 
a solution in TB2 or otherwise the empty set. 

1. cddl {}, eg := Gquo(L/), VG ;= ValG(i/) . 

2. Find the base equations in TB2 by comparing GQ and VG with the correspond- 
ing data in the table. 

(a) if there is no match then return 'Not solvable within the Table'. 

(b) if there is a matching equation Lc, cddl := cddl U {Lc}- 

3. For each Lc G cddl, compute candidate values for the parameters using GQ and 



4. Construct a set cdd2 by substituting parameters by the values determined in 



5. For each Lp e cdd2 check if there exists a GT- Transformation from Lp to Lj. 

(a) if there is a GT-transformation then apply GT to the known solution of Lc 
and return the solution. 

(b) if there is no GT-transformation found return 'Not solvable within the Ta- 



Gquo(L) = {a, a, b, b, c, c}, ValG(L) {(0, 4)} 




VG. 



StepO] 



ble'. 
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6 Improvement 



A similar approach can be applied to higher order operators that are solvable in terms 
of order two. Suppose L4 is a fourth order operator that is solvable in terms of order 
two in D. Then L4 is equal or gauge equivalent to either Kf^ or Ki(§)K2 for some 
second order operators Ki, K2 E D with nonvanishing coefficients. The candidates 
for Ki can be detected analogously using Theorem l3.15l 

Concerning the applications to proving positivity of special functions inequalities 
it has to be noted that representations in terms of finite linear combination of squares 
with non-negative coefficients need not exist on the full range of validity of a given 
inequality, as can be seen below. Further investigations of the applicability of this 
approach as well as an implementation of the above mentioned extension to higher 
order recurrences are ongoing work. 

7 Applications 

Our main motivation to extend finding closed form solutions of difference equations in 
terms of symmetric products is to develop an algorithmic approach for proving special 
functions inequalities. Existing symbolic methods IIT2I ITSl [161 are based on using 
Cylindrical Algebraic Decomposition (CAD) which in several examples has proven to 
be an effective way for proving positivity of sequences that are given only in terms of 
their defining sequences. However, it is sometimes unsatisfiable to have a proof that 
only comes with "True" without any certificate. Some classical proofs of inequalities 
are using rewriting of the given expression as linear combination of easy to verify 
positive objects such as sums of squares. The present work tries to make this approach 
algorithmic. Certainly it will not provide answers for any special functions inequality, 
but it is a first step in a new direction of automatic inequality proving. Below we give 
two examples, one for each of the cases distinguished above, of classical problems that 
can be solved fully or at least partially using the presented algorithm. Note that all 
of these identities stated can be proven easily using existing algorithms for symbolic 
summation. The novelty is the automatic discovery of certain closed form expressions 
for sequences that are given only in terms of their defining recurrence relation. In this 
sense it is comparable to the above mentioned algorithms based on CAD. 

7.1 Clausen's formula 

Proofs of special function inequalities often depend on a variety of classical techniques 
such as argument transformations, integral representations of hypergeometric series 
and many more. For instance in the proof of the Askey-Gasper inequality |i2J, which 
played a key role in the proof of the Bieberbach conjecture by de Branges ||9l, Clausen's 
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formula 



3P2 



-x,x + a + l, ^ 



a+2 
2 



(6) 



— 2F1 



\x + \{a + 1) 



- + 1 

2^ ^ 



entered at a central point. Zeilberger [TO] has shown how this identity can be proven 
using symbolic summation. By means of the algorithm presented here, Clausen's for- 
mula can be discovered entirely automatic. 

The hypergeometric function in (|6) satisfies a third order recurrence that is given 
by the operator L-j, and is too large to be displayed here. It can however be found easily 
common symbolic summation algorithms 1301 [8l [iTl . This difference operator is the 
input for our procedure and we start by determining the local data given by 



Gquo(i3) 
ValG(L3) 



{(0,2) 



2z- 1 ± 2 \Jz'^ - z 
a, 2)}. 



,-20 + l±2 



A table look-up shows that this local data is compatible with [T] in Section 1431 Com- 
paring local data and solving the system mod Z the following candidates for a, 6 and c 
can be found: 

ae{0,i}, 6e{ia, ia+i}, c e {^a + 1, ia + §}. 

There is no term transformation for these operators and an application of HOM shows 
that we obtain a constant map if a = 0, 6 = + ^ and c = + 1. 



7.2 Turan inequality for Hermite polynomials 

The positivity of Turan determinants has been proven for many different families of or- 
thogonal polynomials. The first Turan inequality was formulated for Legendre polyno- 
mials II25I and Szego ll23l has given four different proofs of this inequality. Szwarc ll24l 
has provided a more general approach for proving Turan type inequalities based on the 
mere knowledge of the recurrence coefficients satisfied by the given sequence. Ger- 
hold and Kauers ifTSl have proven and improved this type of inequalities using their 
CAD-based method. The approach presented here does not give a full proof for Turan 
type inequalities, however it gives a representation of the given determinant in sums 
of squares derived from the third order annihilating operator of the determinant. In 
the case of Hermite polynomials this yields a representation that gives positivity in the 
limit for n tending to infinity. 

Turan's inequaUty for Hermite polynomials Hx{z) reads as follows: 

A,{z) = - H,{z)H,+2{z) > 0, n > 0, z e M. 
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Then an annihilating operator of l^x{z) is L/,. := T^ + (2a; + 2 — 4z^)t^ — 4(.t + 2)(.t — 
2z^ + 4)r — 8(1 + x){x + 2)^ and the local data of this operator is 

Gquo(i,i) = {-1 ± ^J-2zH^ + (z^ ± l) 1 ± 2 ^J-2zH^ - 4 z^tj , 
ValG(i^) = {(0,2)}. 

are elements in Gquo(i/i) and these are equivalent 
to —1 ± 2 z^t^ + z^t under ^2, see Definition l3.12l for ^2- Thus the local data of 
Lh correspond to those of the third entry of the table given in Section l4n 

Using the algorithm described above a gauge transformation can be found that ap- 
plied to Ha:{z)'^ yields the following equivalent formulation 

A,(z) - ii/J+i(z) + 2(a; + 1 - z^)H,izf + 2x^HI_,{x). 

This representation gives the positivity of Turan's inequality for 

z e [-Vx + 1, Va; + 1 ], a; > 0. 
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